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We investigate the role of hyperfine mixing in the electromagnetic decay of ground state doubly 
heavy be baryons. As in the case of a previous calculation on 6 — ^ c semileptonic decays of doubly 
heavy baryons, we find large corrections to the electromagnetic decay widths due to this mixing. 
Contrary to the weak case just mentioned, we find here that one can not use electromagnetic width 
relations obtained in the infinite heavy quark mass limit to experimentally extract information on 
the admixtures in a model independent way. 
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Baryon 
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Sh 
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Quark content 


Sh 


J" 


M [MeV] 




{b c}l 


1+ 


3/2+ 


6996 




{b c} s 


1+ 


3/2+ 


7075 




[b c] 1 


0+ 


1/2+ 


6958 




[b c] s 


0+ 


1/2+ 


7038 




{b c}l 


1+ 


1/2+ 


6928 




{b c} s 


1+ 


1/2+ 


7013 



TABLE I. Quantum numbers and quark content for unmixed ground state doubly heavy be baryons with weU defined Sh (spin 
of the heavy quark subsystem). The masses were obtained in Ref. T7]. 5*^ stands for the spin and parity of the heavy quark 
subsystem, while stands for the total spin and parity of the baryon. For J — 1/2, actual physical states are mixtures of the 
Si,c, 5j,^ {^bc, ^bc) states, and they appear in Table HIl 



I. INTRODUCTION 



In the infinite heavy quark mass limit, and according to heavy quark spin symmetry (HQSS) [l|, one can select 
the heavy quark subsystem of a doubly heavy baryon to have a well defined total spin Sh = 0, 1. This has been the 
default assumption by most calculations of doubly heavy baryon spectroscopy f2l-[l^. However, due to the finite value 
of the heavy quark masses, the hyperfine interaction between the light quark and any of the heavy quarks can admix 
both Sh — and 5/i = 1 spin components into the wave function. For ground state (total orbital angular momentum 
L = 0) be baryons, one should expect the actual physical 5 particles (quark content bcu or bed) to be mixtures of the 
Sfec {Sh = 1) and Sj^^ [Sh — 0) states. Similarly in the strange sector, the physical £7 particles (quark content bcs) will 
be mixtures of fibc {Sh — 1) and {Sh — 0)) states. 

While mixing effects are negligible in the spectrum, it was pointed out in Ref. that hyperfine mixing could 
greatly affect the decay widths of doubly heavy baryons. The calculation for & — )■ c semileptonic decay of doubly heavy 
baryons was conducted by the same authors in Ref. [l6l |. where they found that hyperfine mixing in the be states 
had a tremendous impact on the decay widths. We qualitatively confirmed their results in Ref. although our 
predictions for the decay widths were roughly a factor of two larger. There, we also showed how HQSS predictions for 
& — >■ c semileptonic decay, could be used to experimentally obtain information on the admixtures of the be baryons in 
a model independent manner. Unfortunately those ratios involved weak decays that have competing electromagnetic 
(e.m.) decays and thus, they will be difficult to observe experimentally. In this context, it was clear the possible 
relevance of hyperfine mixing effects in e.m. decays. In fact, the authors of Ref. 16] expected hyperfine mixing 
effects to play an important role also for e.m. transitions. As a result of these considerations we included in Ref. [1^ 
predictions for ratios that involved e.m. decay widths evaluated in the infinite heavy quark mass limit. This limit 
implies that the spin of the heavy quark subsystem can not change in an e.m. transition. 

In this letter we perform the full calculation using the same quark model as in Ref. [l3]- To our knowledge there 
is only one prior calculation of e.m. decays of doubly heavy be baryons [Tsj . There, the authors used the e.m. 
radiation as a means to investigate the diquark structure but no hyperfine mixing was considered. In this work we 
restrict ourselves to transitions involving ground state (L = 0) be baryons and our emphasis is put on the relevance 
of hyperfine mixing for those transitions. In Table |T] we show the quantum numbers for the ground state of unmixed 
be baryons classified so that Sh is well defined. The physical spin-1/2 be states are mixtures of the Sbc, SJ,^ (^ibc, ^^bc) 
states shown in that table. Their quantum numbers and admixture coefficients appear in Table |TT1 As for the weak 
& — >■ c decays analyzed in Ref. (l7j . we find here that hyperfine mixing largely affects the e.m. decay widths. On the 
other hand we find that contributions that change the spin of the heavy quark subsystem are very important in the 
evaluation of the e.m. decay widths. We are thus far from the infinite heavy quark mass limit according to which the 
spin of the heavy quark subsystem can not change in an e.m. transition. Due to this fact the e.m. decay width ratios 
proposed in Ref. [l7[, and obtained within that assumption, are not valid for the actual heavy quark masses and can 
not be used to experimentally extract information on the admixtures in a model independent way. 

The paper is organized as follows: in Sec. |ll] we collect general formulas to evaluate the e.m. decay width. We 
also give an appropriate form factor decomposition of the electromagnetic current matrix elements as well as showing 
how the form factors can be obtained in terms of those matrix elements. In Sec. IIIII we present our nonrelativistic 
states and the way the matrix elements are evaluated in our model. Finally in Sec. IIVI we present the results and the 
conclusions of our work. 
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r 
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wW - 

^bc ^ 


0.902 4^+0.431 Ebc 


1/2+ 


6967 




0.899 n'tc 


+ 0.437 Qtc 


1/2+ 


7046 


=■(2) _ 
"6c ^ 


-0.431 H'i,^ + 0.902 Etc 


1/2+ 


6919 




-0.437 n'tc 


+ 0.899 ntc 


1/2+ 


7005 



TABLE II. Physical spin- 1/2 doubly heavy be baryons. The admixture coefficients and the physical masses were obtained in 
Ref. [H. 



II. ELECTROMAGNETIC DECAY 

The electromagnetic decay width for the i? — > B' 7 process is given bjQ 



2M J {2tt)^2E' J (27r)3 2w 



2 J 



s s' r 



where P = (Af;0), P' = (£" = \/ Af'^ + P''^^P') are respectively the four momenta of the initial and final baryons. 
J is the total spin of the initial baryon and s, s' are the spin third components of the initial and final baryons. 
q = {lj = \q\^q) is the final photon four momenta, being £r{q) its polarization vector. Finally .jf^ '^\P,P') stands 
for the electromagnetic current matrix element 



with 



J^s^ '{P,P') = (b', s'P' |J^„(0) \B,sP^Q 



in 



(3) 



(4) 



where the different are the quark charges in units of the proton charge e, and aem is the fine-structure constant. 
Due to the conservation of the electromagnetic current we can take for real photons 



-5" 



(5) 



and thus rewrite the total width as 

1 f d^P' 

r = 



2M J {2ttY2E' J (27r)3 2w 



{2t:)H^^\P - P' - q) 



2 J 



^E T.^^s'^' '^p^p'){^^s?',{p,p') 



(6) 



The double sum in Eq. dH]) is a Lorentz scalar and it can only depend on — IvP , P'^ — M'"^ and P ■ P' = ME' , 
with E' — (M^ + M'^)/2M. As a result, all integrals can be done explicitly and we have the final expression 



1 M^-M'^ -1 



lE H^^s"^' '(P^P')[^^s?',{P,P') 



SttA'P 2M 2 J 

s s' 

where, for the purpose of evaluation, we shall take q along the positive Z— axis. 



(7) 



A. Form factor decomposition of the electromagnetic current matrix elements 

We will analyze 1/2 1/2 and 3/2 1/2 transitions. 



^ Note the normalization of the baryon states should be such that 

(b, s'P' \B, sp) = {2TTf 2E S'-^'' {P - P') 



(1) 
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1. Case 1/2 ^1/2 

For 1/2 — > 1/2 transitions we can write the following form factor decomposition 
^(^, = {b\ s' P' = -q\J!f,^{0) \B,sP = 6 



2(M-M')P'^\ fP^" {M^ - M'^ + q^)P't' 

7 - ITPi 7775 + [ — 7T7T7^ 7T75 — I ^2 



M2 - M'2 - g2 



7\f M(M2 - M'2 - g2) 



Us(0) 



(8) 



where u(0), u'{—q) are the Dirac spinors (normalized to twice the fermion mass) for the initial and final baryon and 
Fi, F2 are form factors that could only depend on the baryon masses and q^. The above form factor decomposition 
trivially satisfies q^J^^^ — 0. For the present case we need the value of the form factors at q2 — q ( \q \ = ^ )• 
We shall have for the double sum in Eq. ([T]) 



2J + 



l9l = ^ 

2 'W M + M' 



M 



M + M 
2P'^ 



/ 2P' 

= 2(M-M')' Pi'|,.,_ m2_m'2 

where in the second equality we have used current conservation. 
The Fi can be obtained as 



m-- 



Fi = 



19 



1 lE' + AP r,^,, 

>^-l/2 1/2l^i^ ) 



2M 



where we have taken q along the positive Z— axis. 



M'2 



M'2 



(9) 



(10) 



2. Case 3/2 ^1/2 

For this case we could use 

Jf/ ^(P,P') = (P', s'P' = -q-\JUO) \B, sP = 6)= u'A-q)f'^''u^ .(0) (11) 



where Ua{0) is a Rarita-Schwinger spinor for the initial spin 3/2 baryon and F"^ is given by 

/ (^V fV fiV 

P"'' -[-^ (.9"^ i - I'^r) + YTa (.9"'' 9 • ^ - 9"^^) + xfe (.9"^ q • P' - g-P"^) ) 75 (12) 



C^, C'l', are vector form factors that, as before, could only depend on baryon masses and q^ . For the double sum 
in Eq. ([7]) We have in this case 



2J + 



^ E E '(^' p') {j^s'^'^iP, P') 



-ipr {(/" + M')f "^(-1)(/' + M)G„^ 7°(r''^)^ 7°} 



19 1 Sin — 



(13) 
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with 



- - -i^ip - 3 ^ + 3 M 

Taking q along the positive Z— axis, the , CY , CY form factors can be obtained as 



(14) 



cY 



cY 



M' 



1 



q\ V 2M(i?' + Af ') ( V2 ^-^=^^2 ^'^ + V 2 ^-''i^^ V2 



Af'^ E' + M' 



M\q\ 



2M 



3 A//^;' - Af'2 



2 Af - 



-^M'J%Y,iP,P')^ 



Af'2 lE' + M 

13 



(15) 



Within our model we shall obtain (see next section) 
so that 



(16) 



cY = G , = , cY 



3 1 



2A/'2 



M ' V 2 |g| y M{E' + M') 

For that case, the trace in Eq. ([T3|) can be evaluated to be 



7T ^-i/2i/2(-F'j-F'') 



(17) 



(Af - M'f [M + Mr ,^v.2 I 



including the —1/4 factor. 



(18) 



III. NONRELATIVISTIC STATES AND MATRIX ELEMENTS EVALUATION 

In this section we briefly describe our nonrelativistic states and the calculation of the electromagnetic current matrix 
elements within our model. 



A. Nonrelativistic states 



Our nonrelativistic states are constructed as 



2Ef 



a2 P2 ^ -=-P + Q2 



asP3^^P-Qi-Q2 ) (19) 



where aj represents the spin (s), flavor (f) and color (c) quantum numbers { a = (s,/, c)) of the j-th quark, and 
{Ef^ , pj), mf- are its four-momenta and mass. M is given by M — ruf-^ + rrif^ + mj^. 
Quark states are normalized such that 



{a' p'\ap) =S^,a, (2^)3 2Ef J^^) (p'-p) 



(20) 



In the transitions under study the baryons involved have bcl quark content, where / represents a light quark u,d, s. 
We choose the wave functions such that quark 1 is a 6, quark 2 is a c and quark 3 is the light one I. 
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^ifo^oaCQi'Qa ) is the internal wave function in momentum space, being Qi (Q2) the conjugate momenta to the 
relative position ri (r2) between the light quark and the b (c) quark. This wave function is normalized as 

Jd^Q,Jd'Q2 J2 (V'ifil3(^i'^2))*V^if.1a3(Qi,Q2) = <^s'. (21) 

ai,a2,a3 

Thus, for our nonrelativistic baryon states we get 

^^(^B,s' P'\B,sP'j^^^5s's (27r)3 S^^^P' - P ) (22) 
For unmixed states with a well defined Sh value the "ipifa^asi Q11Q2) wave function has the general form 

i'ifk'IcAQi^Q^ ) ?''( Qi, Q2 ) Sf, b Sf, c % Kl/2, 1/2, Sh; s,,S2, si + S2){Sh, 1/2, J; s, + S2, S3, s) (23) 

where "^"^J^,"^ is the color wave function, with C3 the fully antisymmetric tensor in three (color) indices, and the 
(ii)i2, j; wi, m2, m) are Clebsch-Gordan coefficients. 

Details on the calculation of the orbital wave function in coordinate space for each of the unmixed states involved 
in this study can be found in Refs. [Tsl. [l9j. 



B. Matrix elements evaluation 



We evaluate the electromagnetic current matrix elements as 

J-f/' ^(^, n = (b', s' P' = -q I Ji;„(0) \B,sP^Q 



em V 

3' „' D' 



'2M^2E' (B\ s' P= -q\J^„AQ) \ B, sP = 0) 



= V2MV2E^ jff ^(g) 



NR 



(24) 



with 



j k 

< ^ (1/2, 1/2, Shj;si,S2,si + S2){Shj, 1/2, J; si + S2,s - si - S2,s) 



(1/2, 1/2, 5,;^; si + s' - s, §2, si + S2 + s' - s){S',,k, 1/2, J'; si + 53 + s' - s, s - si - S2, s') 



e 
3 

2e 



{Qi-^^q,Q2 + =q. 



W 



Ub si+s'-siQi - q) 7'' UbsiiQi) 



2Ebi\Qi - q\)2Eb{\Qi\) 



rrib + mi \ * si+s'-s{Q2 - q) 7^ "csi {Q2) 



M' 



2E,{\Q2 - q\)2E,{\Q2\) 



+eie [4>k\Qi + =q',Q2 + =q')j <5s.,s;,, 

^^(5/,. 1/2, J; rn, s — m, s){Shj, 1/2, j'; m, s' — m^s) 

„ ^' -m(-Qi - Q2- q) 7^ Mi^-m(-Qi - Q2) 1 

2i;/(|-Qi-g2-'f|)2£;;(l-Qi-Q2|) J 



(25) 



where we have used the one-body approximation. The first two terms are the contribution from the b and c quarks 
respectively, whereas the third term is the contribution from the light quark, e; is the charge of the light quark in 
units of the proton charge e. Besides, we sum (sums on j, k) over the different contributions to the physical states 



7 



and the , factors are the corresponding admixture coefficients. For the evaluation of the matrix elements we 
take q along the positive axis. 

In order to be able to evaluate the spin sums explicitly it is useful to use the following relations obtained assuming 
q to be along the positive Z— axis 



V2E'2E 



Us'{p' ^ p - q)j°Us{p) 



{E' + m){E + m) 



2E'2E 



{E' + m){E + m) 



2E'2E 



p2 _ \q\p3 



{E' + m){E + m) 



-X. 



1 



p2 _ |^|p3 



{E' + m){E + m) 



+i 



{E' + m){E + m) 



(E' + m){E + m) 



( {-p^ +ip^)5s's+i +{p^ +ip^)S,,s^i) 



(26) 



where we work in Pauli-Dirac representation and x stands for a Pauli spinor. Similarly for the spatial components 
one has 



^/^uAf = P- q Wu,{p) - sj ^^^^§^XI 



[p-q]' 



+i 



E + m E' + m 
[E' + m)(E + m) [E' + m) J 



{E' + m){E + m) 



2E'2E 



{p- q)^ 



E-E' 



E + m E' + m {E' + m){E + m) 
\q\{E + Tn) - {E-E')p^ 



{-p^Sji +p^5j2){5si/2 - 5s -1/2) ) 5. 



{E' + m){E + m) 



s' s—1 ^s' s+1 



+i6 
+'5,3 



\q\[E + m) - {E-E')p^ 



{E' + m){E + m 
E-E' 



s' s+1 



+ Ss's-l) 



(E' + m){E + rn) 



ii-p^ + ip^)Ss' s+i + (p^ + ip^)Ss' s-i) 



(27) 



Using the above results in Eg. ([25]) it is now easy to see why for 3/2 — > 1/2 transitions we find Ji/2ip.^Pt P') = as 
a result of the orthogonality relations of the Clebsch-Gordan coefficients. Explicit evaluation of the spin sums also 
shows that for 3/2 1/2 transitions J^3^2-i/2(^'^') ~ V^^-?l/2 i/2(^' ^')- This result can be obtained realizing 
that 

^ (1/2, 1/2, 1; si, 32, si + S2)(l, 1/2, 3/2; 

(1/2, 1/2, S*;,; Sl + s' - s, S2, Sl + S2 + s' - s){S',„ 1/2, 1/2; si + S2 + s' - s, s - si - S2, s') Ss' s+i 

= -1= ([(1/2 ® l/2f^ ® l/2],f I I [(1/2 ® 1/2)1 ^ ^j^^z^^ ^28) 

where | [(1/2(8) 1/2)"^ ® l/2]f) is the total spin state of three spin 1/2 particles coupled to total spin J and third 
component s, and a''_^\ is the +1 component of the spin operator of the first particle. Similarly 

^(1, 1/2, 3/2; m, s - m, s){Sl 1/2, 1/2; m, s' - m, s') 5^ s'Js' s+i 

m 

= --^ ([(1/2 ® 1/2)^^. ® l/2]f I a^_^l I [(1/2 ® 1/2)1 ^ ^/sj^/^) (29) 

(3) 

with a_:^l is the +1 component of the spin operator of the third particle. The Wigner-Eckart theorem now gives 
[(1/2 1/2)^^ ® l/2]%\ a^^] I [(1/2® 1/2)1 ^ 1/2]%^^^) 

= V3 ([(1/2 1/2)^^ ® 1/2]'//^ a^f] [(1/2 1/2)^ ® 1/2]!./,%) (30) 
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r(10-* GeV) 


r(10-* GeV) 


' — 'bcu 


' — 'bcu'y 


4.04 




^^'bcl 


3.69 


' — 'hcd ~~ 




4.04 








' — 'hcu " 




105 


^bc ^ 


Q,bc"/ 


20.9 


' — 'bed 


Hftcdl 


50.5 








' — 'bcu ~ 


' ^bcu^ 


0.992 


^'bc ^ 




0.568 


' — 'bed 


■> S6ci;7 


0.992 









TABLE III. Electromagnetic decay widths (in units of 10 * GeV) for unmixed states with a well defined Sh value. 



r(10-* GeV) 






r(10-*GeV) 






r(10-* GeV) 


r(10-* GeV) 


' — 'bcu 


^ "bcu 1 


6.05 


^Ic- 




0.31 


' — 'bcu 


bcu 1 


1.56 


^le- 




0.415 


' — 'bed 


^ "bed 1 


0.12 








' — 'bed 


^ "bed 1 


0.748 








' — 'bcu ' 


^ "bcu 1 


73.9 


^tc- 




50.2 


' — 'bcu ~ 


^ S<^' 7 

beu ' 


123 


^bc ^ 


> ^^rc'7 


24.2 


' — 'bed ~ 

=■{!) 

' 'bcu 

=■(1) _ 

' 'bed 


^(2) 

■> "bcd-y 

^ "bcu 1 
^ "bed 1 


103 
12.4 
20.9 


'''■be 




8.52 


' — 'bed ~ 
' 'bcu 

=■(1) _ 

' 'bed 


^ "bedl 
^ "bed ' 


59.2 
22.8 
11.0 


"be 


^^1^ 


3.78 



TABLE IV. Electromagnetic decay widths (in units of 10~* GeV) for physical states. We show the full calculation results 
(left panel) and results obtained considering only the contributions where the total spin of the heavy quark subsystem does not 
change (right panel). 



IV. RESULTS AND CONCLUSIONS 



In Table Hill we show the results for the e.m. decay widths evaluated with the unmixed states of Table HI while in 
the left panel of Table IIVI the results using the physical spin- 1/2 be states of Table HIl are given. The effects of mixing 
are relevant for all transitions. In particular we find 

^{"bcd ~^ "bcd^) ~ 33 ^{"bcd ^ "bcdl) 

Tini'J ^ ni'Jj) « 15 r(r!;, ^ (31) 

which shows very clearly that hyperfine mixing can not be ignored when evaluating e.m. transitions involving spin-1/2 
be baryons. Besides, we observe that mixing breaks the degeneracy originally present in the unmixed case for most 
transitions involving bcu and bed baryons. 

We know that in the infinite heavy quark mass limit only the terms where the spin of the heavy quark subsystem 
does not change can contribute to the decay widths. This effect can already be seen in Table IIIII in the fact that 
T{El^i — > Ebcn) ^ r('^6c/ ~^ "bcH) '^^ ™ smallness of r(Sj^; — Sbc/7)- Similar results are observed in the V, 
sector. 

To see how far from the ideal infinite heavy quark mass limit we are in this case, we show in the right panel of 
Table HV] the results obtained with the physical spin-1/2 be states, but considering only the contributions to the 
decay widths of the terms where the total spin of the heavy quark subsystem does not change. We see big changes 
when compared to the full results in the left panel of the same Table. As a result, and in contrast to the weak 
decay case discussed in Ref. (l7| . and to our prior expectations also outlined in this latter reference, heavy quark 
spin symmetry relations deduced in the infinitely heavy mass limit, are not accurate enough for the study of e.m. 
transitions involving doubly heavy baryons. Next-to-leading corrections turn out to be quite large as the differences 
between the two panels in Table IIVI indicate. Another important point here is that the decay widths are proportional 
to the factor (M^ — M'^)(M — M')^ coming from phase space and the spin sums. As M is close to M', the decay 
widths are very sensitive to the actual baryon masses. 

Finally, we would like to stress, once more, that the experimental measurement of e.m. widths will be extremely 
valuable in order to extract information on the hyperfine mixing of doubly heavy be baryons, as the difference among 
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the results of Table IIIII (for unmixed states) and the left panel of Tabic IIVI (mixed states) clearly show. However, 
we should also point out that by looking only at e.m. transitions, it would not be possible to determine their actual 
mixing matrix without relying on a theoretical model. In this respect, the situation is more favorable in the case of the 
semileptonic weak decays of these baryons, as we discussed in Ref. [17|, where leading order heavy quark symmetry 
relations turned out to be much more accurate. 
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